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A two-component system is considered, which is subject to external and possi¬ 
bly fatal shocks. The lifetimes of both components are characterized by their 
hazard rates. Each shock can cause the immediate failure of one single or of 
both components. Otherwise, the hazard rate of each component is increased 
by a non fatal shock of a random amount, with possible dependence between 
the simultaneous increments of the two failure rates. An explicit formula for 
the joint distribution of the bivariate lifetime of the two components is pro¬ 
vided. Some positive dependence properties of the bivariate lifetime are found 
out. The influence of the shock model parameters on the bivariate lifetime is 
studied. 

Keywords: Bivariate non-homogeneous compound Poisson process, Haz¬ 
ard rate process, Aging properties. Positive dependence properties. 


1. Introduction 

This paper is devoted to the survival analysis of a system subject to compet¬ 
ing failure modes within an external stressing environment. The external 
environment is assumed to stress the system at random and isolated times 
according to a random shock model. Such a model can represent external 
demands e.g., which put some stress on the system at their arrivals. The 
occurrence of shocks is classically modeled through a non-homogeneous 
Poisson process and components lifetimes are characterized by their failure 
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rates. The shocks are simultaneous for both components and each shock 
can be simultaneously fatal to both components, which induces some depen¬ 
dence between the components. Also, a shock increases the failure rates of 
the surviving components of a random increment, with possible dependence 
between simultaneous increments. Finally, following [1], the probability for 
a shock to be fatal depends on the shock’s arrival time, which induces a 
last kind of dependence. Based on this setting, the aim of the paper is the 
study of the bivariate lifetime of the two components. 

The paper is organized as follows: the model is specified in Section 2. 
An explicit formula for the joint survival function of the bivariate lifetime 
is provided in Section 3, as well as some positive dependence property. The 
influence of the model parameters on the bivariate lifetime is also discussed. 
Numerical experiments hnally illustrate the study in Section 4. 


2. The model 

Out of the stressing environment, the two components are assumed to be 
independent. The lifetime of each component is characterized by its intrin¬ 
sic hazard rate hi(t), i = 1,2, or by the corresponding cumulative hazard 
rate Hi{t) = Jq h i{u)du, i = 1,2. Stresses due to the external environment 
arrive by shocks, independently of the system intrinsic deterioration. The 
shocks occur at time Ti,T 2 ... according to a non-homogeneous Poisson pro¬ 
cess (fVt)t>o with intensity X{t) and cumulative intensity A(t) = f* X{x)dx. 
The shock at time T„ increases the hazard rates of both components of 
a random amount Vn = {Vn^\ , where the increments are in¬ 

dependent and identically distributed (i.i.d.) and independent of the shocks 
arrival times {Tn)n>i- The simultaneous increments and Vn^'^ at time 
Tn may however be dependent. Furthermore, a shock can be fatal, and can 
possibly induce the immediate failure of one single or of both components. 
The fatality of a shock does not depend on the system intrinsic deteriora¬ 
tion but depends on the shock’s arrival time. The following notations are 
used: 

• Poo (Tn)' probability that the shock at time r„ induces the simul¬ 
taneous failure of both components, 

• Pii (Tn): probability that the shock at time r„ induces no failure 
at all among the two components, 

• Poi (Tn): probability that the shock at time T„ is fatal only for the 
first component. 



• pio {Tn)' probability that the shock at time T„ is fatal only for the 
second component, 


with Eo<*j<iKi (•) = 1- 

The common distribution of the i.i.d. random vectors Vn = , 

n G N* is denoted by ^(dvi,dv 2 )- When subscript n is unnecessary, we 

drop it and set V = to be a generic copy of V„ = . 

For j = 1,2, the distribution of jg denoted by ^j{dvj). 

We set (At)t>o = ( ) to be the bivariate compound Poisson 

“ V /t>o 

process defined by 

Nt / Nt Nt \ 

n—1 \n=l n—1 / 

where En=i ■■■ = 0 and T = (t{As,s > 0) is the cr-field generated by 

Providing that it is functioning up to time t, the conditional hazard rate 
of the {i = 1, 2) component at time t given = a{As, s > 0) is 

Vzit) = hi{t) + 


Let Ti, i = 1,2 be the lifetime of the component without taking into 
account the possibility of fatal shocks and let z = 1, 2 be the time of the 
first fatal shock for the component. We have 

= ( 1 ) 


and 


Nt 

fc=i 

where qi (•) = pn (•) + pio (•) and q 2 (•) = Pn (•) + Poi (•)■ 

Given T, the random variables Ti and T 2 are assumed to be condition¬ 
ally independent one with the other, and conditionally independent of 
and ^ 2 - However, as the two components may fail simultaneously at each 
shock, and ^2 are not conditionally independent given T. Given T, their 



conditional joint survival function is 




( Nt Ns 

n^ii(T*) qi{,Ti) if s>t, 

Ns Nt 

]^Pii(T*) g2(T*)ifs<t 

>. 2=1 i=Ns-\-l 


( 2 ) 


where JlLi - = 1- 

We set Y = (Id) ^ 2 ) to be the lifetime of the two components. It is easy 
to see that Yi = min(Ti, ^i) for i = 1,2. 


3. Theoretical results 

In all the paper, the results are provided without proof, due to the reduced 
size of the paper. 


Proposition 3.1. The joint survival function ofY = (Yi,!^) is given by 
FY{s,t) 


— g-Hi{s)-H2{t)-A(meix{s,t)) 


i(s,t) 


/2((s — w)~^, {t — w)^)pii{w)X{w)dw 
\Jo 

^ h{s - w)qi{w)X{w)dw + il 2 {t - w)q 2 {w)X{w)du?j 

(3) 

where ft stands for the Laplace transform of the bivariate distribution p of 
V, with 

fl{xi,X 2 ) = JJ p{dvi,dv 2 ) for all xi, X 2 > 0 

and stands for the Laplace transform of the distribution of V^'‘\ 
with 


Pi{xi) = J e pi{dvi) = E, ^e for i = 1,2 and all Xi > 0. 

We recall that Y 2 is said right-tail-increasing in Yi, written RTI{Y 2 \Yi), 
as soon as 


P(Y2 > a;2|Yi > a;i) 

is non decreasing in xi for all X 2 > 0. Also, Y 2 is said left-tail-decreasing in 
Yi, written LTD{Y 2 \Yi), as soon as P(Y 2 < a; 2 |Yi < Xi) is non increasing 
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in xi for all X 2 > 0. Both i?T/(F 2 |hi) and LTL)(Y 2 lVi) properties are pos¬ 
itive dependence properties, which imply association and positive quadrant 
dependence of Y, see [5] for more details on these different notions. 

Proposition 3.2. Both properties RTI{Y 2 \Yi) and RTI{Yi\Y 2 ) are true. 

However, as will be shown in the next section, it is possible that one 
margin is not left-tail-decreasing in another margin, that is LTD{Y 2 \Yi) 
property is no always true. 

We now study the influence of different parameters on the bivariate life¬ 
time Y. With that aim, two similar systems are considered, with identical 
parameters except one. We add an upper bar to all quantities referring to 
the second system. (For instance, we use A (w) for the second system). The 
next result shows that, as expected, the more frequent the shocks are, the 
smaller the bivariate lifetime Y is. 

Proposition 3.3. Let us consider two systems with the same parameters 
except from the intensity of the non-homogeneous Poisson process. Assume 
that X{w) < X{w) for all w > 0. Then, Y is smaller than Y in the sense of 
the Weak Hazard Rate order (Y <whr Y), whieh means that FV(x)/.Fy(x) 
is increasing inx.= (xi,X 2 ) € {y : Fy (y) > 0}). As a eonsequence [f], we 
also have Y <uo Y (Upper Orthant Order), namely Fy{^) Y Fy{^) for 
all X e and 

my(x) = E(y — x|F > x) < my(x) = E(y — x|F > x) 

for all X G where my and my are the multivariate mean residual 

lifetimes of Y andY. 

The next result shows that the smaller the fatality of shocks is, the 
larger the bivariate lifetime Y is. 

Proposition 3.4. Let U (w) with distribution P {U (w) = {i,j)) = Pi j (w) 
for all i,j G {0,1} and U (w) defined in the same way with respect to 
the family {pij Let us assume that U (w) <uo U (w) for all 

w > 0. Then, Y <uo Y. 

Finally, assume that V and V are marginally identically distributed. In 
that case, Y and Y also are identically distributed. The upper orthant order 
is then equivalent to the Positive Quadrant Dependence order (PQD), which 
compares the dependence between the margins. We recall that the PQD 
order implies the Laplace transform order [3] so that the Laplace transform 
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Figure 1. Example 4.1, RTI and not LTD property 


order can also be considered as some comparison of the dependence between 
the margins. In that setting, the following result means that the more 
dependent the failure rates increments are (in the sense of the Laplace 
transform order), the larger the bivariate lifetime Y is. 

Proposition 3.5. Let V and V be marginally identically distributed. As¬ 
sume V to be smaller than V in the sense of the bivariate Laplace transform 
order (V <cV), namely E < E 

Si, S 2 > 0. Then Y <pqd Y. 

4. Numerical experiments 

In all the following experiments, we take /li = 0 for z = 1, 2. Also, we take 
pld = l/Yl + z = 1 , 2 , where and are independent. 

Example 4.1. The parameters are: A(a;) = e®,pii(x) = poi{x) = pio{x) = 
Also, the IfYl’s (z = 1,2,3) are exponentially distributed with respec¬ 
tive parameters 1, 5 and 6 . Figure 1 (left) shows the right tail with 

respect to x for various values of X 2 . Whatever X 2 is, we observe that the 
right tail is always increasing so that, as expected, the RTL iY 2 \Yi) prop¬ 
erty is true. However, the left tail been plotted in Figure 

1 (right) and we observe that it is not monotonous. In that example, the 
LTD (I 2 IL 1 ) property is consequently false. 

Example 4.2. The parameters are: Pii(a;) = Poi{x) = Pio{x) = 
and the [/(*^’s (z = 1,2,3) are gamma distributed with parameters (oj, 1) 
and ( 01 , 02 , 03 ) = (1,2,3). We consider A(a;) = x and A(a:) = 2x, with 
A(x) < A (a;). The quotient r(x) = Fy{x)/Fy{^) is plotted in Figure 2, 
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Figure 2. Example 4.2, t(x) is monotone 



Figure 3. Example 4.3, Influence of fatality of shocks 

where V and Y refer to A(a;) and A(a;), respectively. We observe that, as 
expected, r(x) is increasing in x. It means Y <whr Y. 

Example 4.3. This example shows numerically the influence of the fatality 
of shocks on lifetime. Let Y and Y be the lifetimes corresponding to p = 
{P^i = l))^,j=o^ and p = {pn = i,poi = g,.Pio = I,Poo = ^)- The E^^^’s 
are the same as in Example 4.2. The other parameter is A(a;) = 2x. We 
have p <pQD P and as expected, the difference d = Fy — Fy is always 
positive (see Figure 3). It means the bivariate lifetime Y is larger than Y, 
in the sense of the positive quadrant dependence order. 



Figure 4. Example 4.4, Influence of the dependence between y(l) and 1/(2) 

Example 4.4. This example illustrates the influence of the dependence 
between and on the bivariate lifetime. All parameters are similar 
as those in Example 4.2, except from the fact that, for the parameters of 
the gamma distributed i = 1,2,3, we take ( 01 , 02 , 03 ) = (1,2,3) and 
( 01 , 02 , 03 ) = (3,4,1). Then, and (resp. and are iden¬ 

tically gamma distributed with parameter (4,1) (resp. (5,1)). However, it 
is easy to check that V <c V. We observe in Figure 4 that, as expected, 
the difference D = Fy — Fy is always positive, so that Y is larger than Y, 
in the sense of the upper orthant order. 
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